Magnetotransport in lateral superlattices with small-angle impurity scattering: 

Low-field magnetoresistance 
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An analytical study of the low-field magnetoresistance of a two-dimensional electron gas subject to 
a weak periodic modulation is presented. We assume small-angle impurity scattering characteristic 
for high-mobility semiconductor heterostructures. It is shown that the condition for existence of the 
strong low-field magnetoresistance induced by so-called channeled orbits is rf^ 2 ql 3> 1, where rj and 
q are the strength and the wave vector of the modulation, and I is the transport mean free path. 
Under this condition, the magnetoresistance scales as rj 7 ^ 2 . 



I. INTRODUCTION 

The effect of a periodic modulation (lateral superlat- 
tice) on the transport properties of a two-dimensional 
(2D) electron gas has been the subject of intensive re- 
search during the last decade. This interest was trig- 
gered by an experiment of Weiss et alM who discovered 
that a weak one-dimensional modulation (grating) with 
a wave vector q || induces strong commensurabil- 
ity oscillations of the magnetoresistivity p xx (B), while 
showing almost no effect on p yy (B) and p xy (B). The 
oscillation minima have been found to satisfy the con- 
dition 2R c /a = n — 1/4 with integer n, where R c is 
the cyclotron radius and a = 2n/q the grating wave 
length. While these findings can be explained in terms 
of the quantunwmechanical band structure induced by 
the modulationB, the phenomenon is, in fact, of qua- 
siclassical nature, as was demonstrated by Beenakkerp. 
He showed that the geometric resonance of the cyclotron 
motion in the grating induces an extra contribution to 
the drift velocity of the guiding center, whose root-mean- 
square amplitude oscillates as \cos(qR c — 7r/4)|. These 
arguments were corroborated by an analytical solution 
of the Boltzmann equation within an expansion in the 
relative strength r/ of the modulation, the result for the 
modulation-induced oscillatory magnetoresistivity Ap xx 
being pJi-fhe order ij 2 . However, although the theoretical 
resultsM accounted for the above experimental features, 
they disagreed strongly with the experiment as far as the 
damping of the oscillations with decreasing magnetic field 
is concerned. The reason for this was an oversimplified 
treatment of the impurity scattering: while the theory 
ofau assumed isotropic scattering, in experimentally rel- 
evant high-mobility samples the random potential is very 
smooth, so that the scattering is of small-angle charac- 
ter, with the total relaxation rate r s _1 much exceeding 
the momentum relaxation rate r _1 . This gap in the the- 
ory was filled in by our previous papeaj where the effect 
of small-angle scattering on the Weiss oscillations was 
studied analytically. It was found that the small-angle 
scattering changes the dependence of the oscillation am- 
plitude on the magnetic field B completely, leading to a 



much stronger damping of oscillations with decreasing B, 
in very good quantitative agreement with experimental 
data. 

This paper continues our investigation of the effect of 
small-angle scattering on the transport in laterally mod- 
ulated structures. We will address the issue of the low- 
field magnetoresistivity, which has been left aside in our 
earlier papero. In combination, Ref. ^ and the present 
paper provide a complete quasiclassical theory of mag- 
netoresistivity of a 2D electron gas subject to a weak 
one-dimensional modulation and a smooth random po- 
tential. 

A distinct low-field magnetoresistivity was observed, 
along with the cflmmensurability oscillations, in the orig- 
inal experiments, as well as in numerous later experi- 
ments on the transport in a lateral superlattice. Specif- 
ically, in low magnetic fields B a positive magnetoresis- 
tivity was found, followed by a maximum in p xx (B). For 
not too strong modulation, the relevant magnetic fields 
are much weaker than those where the Weiss oscillations 
are observed, so that the two effects can be easily sep- 
arated. Soon after the first experimental observation it 
was understood^ that the low-field magnetoresistivity is 
related to the existence of open (channeled) orbits in the 
magnetic fields B < B c = (r)c/2e)qm,VF. It is worth 
mentioning that this effect ,which is not found within the 
^-expansion used in Refs. [Sflj, has its counterpart in the 
context of the sound absorption in metals in the presence 
of a magnetic field. There, the trapping of electrons in 
channeled orbits by a sound wave leads to non-linearity of 
the acoustic response of an electron gas, as-was observed 
experimental^ and,analyzed theoreticallyQ. Though the 
work by Beton et alB (see also a more recent papem) ex- 
plained qualitatively the low-field magnetoresistivity as 
an effect of channeled orbits, a quantitative analytical 
description of the problem requires that the nature of 
disorder be taken into account. This is done below in 
this paper, where we demonstrate that the small-angle 
character of scattering leads to a parametrically differ- 
ent magnitude of this effect. We will also study how the 
contribution of drifting orbits to Ap xx is modified in the 
low-magnetic-field range, B ~ B c . 
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II. LOW-FIELD MAGNETOTRANSPORT IN 
SMOOTH DISORDER 

A. Generalities 

We consider a weak periodic potential 

V(x) = 7]Epcosqx , ))C 1 (1) 

acting on a 2D electron system (Fermi energy Ef) at 
temperature T = 0. In the presence of a perpendicu- 
lar magnetic field B, the electron at the Fermi level will 
perform cyclotron motion (frequency lu c — eB/mc, cy- 
clotron radius R c — vp/uj c , where vp is the Fermi veloc- 
ity). The periodic potential causes a slow drift motion in 
the y-direction with velocity v,i. The classical equations 
of motion are (the electron charge is — e) 



my 



e. R dV 
—yB- — , 

c dx 

l -xB . 



(2) 
(3) 



Averaging (|2|) over one cyclotron revolution one finds the 
drift velocity 



c UJ C 

eB 7r 



dxdV 
x dx ' 



(4) 



where a; m i n , a; max are the classical turning points of the 
periodic motion. Integration of (||) yields 



y = uj c (x - x Q ) , 



(5) 



where xq is the point on the trajectory where y = (i.e. 
the guiding center coordinate). The a;-componcnt of ve- 
locity is now easily found from energy conservation, 



-(E-V eS (x)) 



1/2 



(G) 



where the effective potential of the motion in x-direction 
is given by 



V cB (x) = V(x) + -^(i - x ) 2 



(7) 



For weak magnetic fields, qR c <C 1, the effective poten- 
tial is given by a parabola centered at xq, with rapid 
periodic oscillations superposed (Fig. |l]a). Intersection 
of this curve with a line of constant energy yields the 
turning points of the classical motion. There is always 
a central allowed region x m - m < x < x max , for which 
the motion consists of a drift of complete cyclotron or- 
bits. For sufficiently strong periodic potential, there ex- 
ists also additional classically allowed region, to the left 
(right) of the turning point a; m i n (x max ). These are so- 
called channeled orbits, winding in a snake-like fashion 
along y. Using the fact that x — xq ~ R c in the region of 



channeled orbits, one finds the condition for the existence 
of channeled orbits B < B c , where 



e 7]mvFq 



(8) 



We will assume that ui c t 3> 1 in the relevant range of 
magnetic fields, which implies for B ~ B c that 77 3> 2/ql, 
where I — vft is the transport mean free path. Since 
ql ~ 300 -T- 1000 in a typical experiment with a high mo- 
bility sample, this condition is fulfilled even for a very 
weak (77 of order of few percent) modulation. 
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FIG. 1. (a) Effective potential V c «{x), Eq. (Rh (b) effective 
potential linearized near a turning point, Eq. (|15|). 



B. Contribution of drifting orbits 

We will first discuss the contribution of the drifting cy- 
clotron orbits to the transport. Substituting (^), (]?]) into 
(I), we find 



Vd = tt v f 

Z7T 



du sin 7 



I — 77 cos u — 



(U - Uq) 1 



21 -1/2 



(qRc) 2 



(9) 



where u = qx, uq — qxo, etc. In the limit r\qR c <C I 
(which is equivalent to B 3> B c ) the term 77 cos u in (||) 
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may be neglected. The integral may then be done, yield- 
ing 



Vd = -T]v F qR c J (qR c ) sinu , 



(10) 



where Jq{x) is the Bessel function. The drift velocity is 

oscillating with uq, i.e. the initial conditions. In order 

to determine the modulation-induced correction to the 

diffusion coefficient, one has to evaluate the average of 
0,2 



2\ — 



2 7 2, 



2^ jo 

—rfv 2 F qR c cos 2 (qR c - tt/4) 

47T 



(11) 



(in the second line we used the condition qR c ^> 1). n Thc 
result ( |ll|) obtained for the first time by BeenakkerB in- 
duces an oscillating (with qR c ) correction 5p xx to the 
resistivity along the modulation wave vector (Weiss oscil- 
lations) . The amplitude of Weiss oscillations for the case 
of white-noise disorder was calculated inH and for the (ex- 
perimentally relevant) case of a smooth random potential 
iro. Here we are interested in the behavior of Ap xx in 
the region of sufficiently low magnetic fields, where the 
oscillations are exponentially damped. This means that 
the position uq of the guiding center changes due to im- 
purity scattering by an amount Suq 3> q^ 1 within one 
cyclotron revolution. Taking into account that (in view 
of the oscillatory factor sin u) the main contribution to 
the integral (J^j comes from the regions close to the turn- 
ing points u m in, Umax, one easily realizes that an average 
squared drift velocity in one cyclotron revolution is then 
given by non-oscillatory part of Eq. (pT|), 



r) 2 v F qR c 



8tt 



(12) 



The diffusive process in y-direction with the velocity 



squared 



and the time step At = 2tt/uj c induces a 



correction to the diffusion constant 



&D yy = (vi) — 

U! r . 



and thus the resistivity correction 



Ap x 



po 



AL> yy = -„ 2 i 



(13) 



(14) 



where po is the Drude resistivity in the absence of modu- 
lation. Eq. ([l4]) reproduces the result of Refs. HQ for the 
saturation value of Ap xx /po at low magnetic fields. 

Having demonstrated how the result ( |l4] ) of 77- 
expansion of the Boltzmann equation is reproduced 
within the present approach for B B c , we are pre- 
pared to turn to the question of our main interest here, 
i.e. to the range of lower magnetic fields, B ~ B c . Since 
the main contribution to the drift velocity, Eq. (^J) , comes 



from the vicinity of the turning points, one can linearize 
the second (parabolic) term in the effective potential (Q) 
near the points x — x$± = xq ± R c , 



V e s(x) - E F = V(x) ± mv F uj c (x - x Q ±) 
= r/E F [cosu ± /3(u — Uo±)] 

We have introduced here the parameter 



(15) 



(16) 



The contribution of the lower limit takes then the form 

1/2 



mm 
v d 



r]v F ( qR c , 

— — ) F(u o ,0), 



(17) 



where 



F(uq, 0) = J du sin it ( u - //,, — — cc» u 



-1/2 



(18) 



The upper limit contribution is found to be ^™ ax — 
— v d nu \u ^-u - The average value of the total drift ve- 
locity Vd = w™ ax + vf™ is therefore equal zero, (vd) = 0, 
where, as before the averaging goes over the position of 
the guiding center, 



2tt jq 



dug . 



The average of the square of Vd is obtained as 

(v 2 ) = 2[(KT ) 2 > - (vT n ) 2 } , 



where we used ((v% 



min \ 2 \ 



((v n d 



max\2\ 



and (v 



mxn ,max\ 



(u™ ln )(wjf ax ). The latter relation is ensured by the impu- 
rity scattering leading to a large shift Suq q^ 1 within 
one cyclotron revolution, as has been discussed above. 

Using (|l^), (pT[), we obtain thus the following contri- 
bution of the drifting orbits to the resistivity correction 



PO 



2ttu! c t- 



2tt 



qlF nc (B/B c ) 



(19) 



where the dimensionless function F nc ({3) is defined by 

F nc ((3) = (F 2 (u , /?)) - (F(u , P)) 2 . (20) 

The function F nc {f3) is shown in Fig. || . In the limit 
(3 — ► 00 one easily finds i ? nc (/3) — > tt/2, reproducing 
Eq. P). 
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FIG. 2. Function F nc (f3) characterizing the magnetic field 
dependence of the contribution of non-channeled orbits to 
resistivity. The dashed line indicates the asymptotic value 
F nc (/3 ^> 1) = 7r/2 corresponding to the result ( |l4| ) of the 
jy-expansion. 



C. Contribution of channeled orbits 

We now turn to the contribution of channeled orbits. A 
particle will spend an average time r c h in a channeled or- 
bit and during this time will propagate with velocity vf- 
If the fraction of particles in channeled orbits is P c h(/3), 
the contribution to the diffusion coefficient will be 



(21) 



In order to calculate P c h(/3) (see e.g. Refs. ||,|8|), we con- 
sider again the effective potential V e e(x) linearized near 
a turning point, Eq. (15). The fraction of the phase space 



(x, 4>) occupied by the channeled orbits is clearly given 
by 



Xc dx 24> max (x) 



(22) 



where x m is the position of a maximum of V e s, x c is the 
nearest point to x rn satisfying V g(x c ) = V e s{x m ) (see 
Fig. |l]b), and 4>max(x) is the limiting angle of the veloc- 
ity vector with the y-axis. In view of rj <C 1 we have 
0max < 1, so that 



Vrn - Vcffjx) 

mvp/2 



1/2 



^max(^) 



Substitution of (E3) into (B2h yields 



V m = V e s(x m ) 



P ch (f3) = -(2r ? ) 1 / 2 $(/?), 



(23) 



(24) 



$(/3) = 



arcsin /3+27T 



du 



4 a/2 J arcsin (3 

Re[y/1 - (3 2 + [3 arcsin (3 - cos u - (3u] 1/2 . (25) 



The function $(/?) satisfies $(0) = 1 and $(1) = and 
is shown in Fig. |3|a. 

The life time of a particle in a channeled orbit, r c h, 
may be obtained by considering the diffusion in the 
space (x, x) in the potential well formed by the potential 
V e f[(x). Particles diffusing beyond the borderline marked 
by the maximum of V c g(x) escape the confining well. The 
lifetime of particles may be defined from the time depen- 
dence of the total number of particles in the well, n w (t), 
as 



Tch 



dtn w (t) , 



(26) 



where it is assumed that at t = there is one particle in 
the well, n w (0) — 1, spread uniformly over the volume f2 
of the corresponding phase space. In order to calculate 
n w we define the phase space density /(x, x, t), subject to 
the initial condition /(x,i,0) = 1/51 and to the bound- 
ary condition / = at the boundary of the confinement 
region. 

To determine f(x,x,t) analytically, we approximate 
the potential well by a parabola 



V eS (x) = AAE 



[x - (x c + x m )/2] 



const 



(27) 



Xrn ^ X *C X c . 



where AE is the depth of the well. Introducing dimen- 
sionless variables 



X 



2[x - (x c + x m )/2] 



Y 



1/2 



and e = E/ AE, the energy takes the isotropic form 



x, (28) 



X 2 + Y 2 , 



e < 1 



(29) 



The scaled phase space trajectories are thus circles. Let 
us now consider the effect of impurity scattering. The 
small-angle scattering by the smooth random potential 
induces a diffusive motion in the space of the velocity 
angle, ((5<fi) 2 ) — 2D$8t, with the diffusion coefficient 
Drf, = 1/t, where r is the transport time. For chan- 
neled orbits we have x ~ 4>vf, implying a diffusion in 
x with diffusion coefficient D± = v f /t. Provided the 
motion of the particles in the potential well is rapid as 
compared to the escape process (requiring the condition 
rfl 2 ql 3> 1, see below), the effect of the anisotropic dif- 
fusion (along Y) is equivalent to isotropic diffusion with 
diffusion coefficient 



2 AAE t 



(30) 



where 



We are thus left with the diffusion equation in a circle, 



4 



DV 2 n(R,t) ~ ^ = , 0< R= (X 2 +Y 2 ) 1/2 < 1 , 

(31) 



supplied with the boundary condition n(l,t) = and 
the initial condition n(R, 0) = 1/tt. The solution is eas- 
ily found by expanding n(R, t) in cigcnfunctions Jo(K n R) 
of the Laplace operator, 



9 i 

7T ^ ^ 



71" Z - ; K n Ji(K n ) 



J {K n R)e 



-Dntt 



(32) 



where k„ are zeros of the Bessel function Jq(x). We 
find thus for the integrated staying probability n w (t) = 

2tt J^dRRn(R,t) 



DKZt 



(33) 



The lifetime of a particle in the channel is now found 
from (PQ) as 



where 



Tch = — V \ = — = -rrr^x (ft) , (34) 



$i(/3) = y/l -fi 2 +/3arcsin/3- |/3 . (35) 
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FIG. 3. a) Functions <J>(/3) and $i(/3) characterizing the 
fraction of open orbits and the life time in a channel; b) func- 
tion _F c h(/3) describing magnetic field dependence of the con- 
tribution of channeled orbits to the resistivity, Eq. (pa). 



The dimensionless function &i(ft) decreasing monotoni- 
cally from 1 at ft = to at ft = 1 is shown in Fig. ||a. 

Collecting the results for P c h(-B) and r c h we find the 
correction to the resistivity due to transport in channeled 
trajectories as 



Ap ch 

rxx 

Po 



2(uj c T) 2 P ch (ft)T ch /T 

^(ql) 2 F ch (ft) , 



where 



F ch {ft) = ft 2 ® x {ft)${ft) 



(36) 



(37) 



Therefore, the magnitude of the low-field magnetoresis- 
tance (i.e. the value of Ap^/po at maximum) scales 
as rf' 7 1 2 with the modulation amplitude. The function 
-Fch(/3) approaches zero for ft — » and ft — > 1 and has a 
maximum of ~ 0.03 at ft ~ 0.4 (see Fig. ||b). 



D. Resistivity correction at zero magnetic field 

Formulas ([l9|) , (|36|) give a zero resistivity correction in 
the limit ft —* 0. However, these formulas are not valid 
for too small ft, and in particular, at B = 0. This is 
clear already from the fact that we used the condition 
lj c t S> 1 in the course of their derivation, which im- 
plies ft ^> 1/rjql. (In fact, the lower limit for ft is even 
more restrictive, see Eq. ([39]) below). Therefore, the zero 
magnetic field case requires a separate treatment. Our 
approach, based on consideration of a disorder-induced 
diffusion in the space of electron orbits, can be, however, 
generalized to the case B = as well. The calculation 
presented in Appendix yields the following result: 



&Pxx{B = 0) = 1/2 



Po 



(38) 



where C ~ 1.24 is a numerical coefficient. How does 
Eq. ([38]) m atch the. finite- B result (|36|)? In our deriva- 
tion in Sec. [IB, QC we assumed that the cyclotron mo- 
tion away from the domain of the phase space occupied 
by the channeled orbits is faster than the escape process 
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from a channel. Since a characteristic angle 
the channeled orbit is 
i 1 / 2 



for 



fy 1 / 2 , we get the condition 



T c h ^ y' /wc, which can be rewritten as 



(39) 



Note that the r.h.s. of (B9J) is assumed to be much 
smaller than unity, see Sec7|IE|. In the opposite case, 
f3 <C l/rj 3 / 2 ql, the zero- field result ( J38| ) should be valid. 
Indeed, comparing ( |36|) and (|38|), we see that they match 
at /3~ l/rfl 2 ql. 



E. Conditions of applicability 

Let us discuss conditions of validity of our consider- 
ation. First of all, the picture of the diffusion in the 
velocity space used above is justified if the particle un- 
dergoes many scattering events within the time r c h ~ ryr, 
which implies 



f] > t s /t. 



(40) 



Assuming a value r/r s ~ 50 characteristic for high mobil- 
ity samples, this inequality is reasonably fulfilled for mod- 
ulation strengths rj > 5%. In the opposite case 77 <C t s /t 
(i.e for very weak modulations) one should replace r c h by 
t s in Eq. (36), yielding 



a ch 2 3 / 2 B 



(41) 



where Z s = %r s , so that the rj 7 ^ 2 behavior changes into 
the jfl 2 one. As to the zero-B formula we expect 
that it transforms at 77 <C t s /t into Ap/p^.^ rfl 2 T/T s , 
matching the result Ap/po ~ 77 s / 2 obtainecB in the limit 
of isotropic scattering (r s = r). 

Another essential assumption was that the oscillatory 
motion within the potential well takes place on a time 
scale much shorter than the escape time T c h. In other 
words, we assume that the particle performs many oscil- 
lations in the channel (with the frequency aj c h ~ rf~^VF<j) 
before it escapes into a conventional cyclotron orbit. The 
condition for this is w c hT c h S> 1, which means that 



V 



3/2 ql > 1. 



(42) 



For experimentally relevant values of ql this implies 77 3> 
1 4- 2%, which is fulfilled in the majority of experiments. 
In the opposite case, t] 3 ' 2 ql -C 1, disorder scattering dom- 
inates over modulation-induced effects, so that a particle 
escapes from the channel before it "recognizes" that it 
is trapped there. In this limit the non-perturbative in 77 
effects related to existence of channeled orbits become ir- 
relevant, and the magnetoresistivity is given by the result 
of the perturbative expansion (i.e. by Eq. ( |l4| ) in the re- 
gion of exponentially damped oscillations), without any 
pronounced features around B = B c . 



Note that in the assumed regime ( fl2| ) the condition 
^ch^ch 3> 1 is still violated in a narrow vicinity of /3 = 1, 
namely at 1 — f3 < (rf/ 2 ql)~ A / 7 , since both w c \i and r c h 
vanish when f3 — > 1, r c h oc (1 — /3) 3 / 2 and w c h cx (1 — /3) 1 / 4 . 
This implies that channeled orbits get gradually de- 
stroyed by disorder as (3 approaches unity and leads to a 
smearing of singularity in Ap xx /po at /3 = 1 [Eqs. jis|), 
(p6|)] over a narrow interval 8/3 ~ {rf'/ 2 ql)~ i / 7 . 

Finally, we discuss the overall magnitude of the ef- 
fect. At B > B c only the drifting orbits exist, and 
Eq. ( fl9| ) predicts an enhancement of Ap xx /p [with re- 
spect to its value well above B c , Eq. (|i~4|)] by the fac- 
tor F nc (/3)~ 1. At B < B c comparison of Eq. @) 
and Eq. (H9) shows that the contribution of channeled 
orbits dominates, yielding an enhancement of Ap xx /po 
[again with respect to (jj"4|)] by a parametrically large fac- 
tor (2 5 / 2 /7r 2 )F(/3)r7 3/2 g/. Let us note, however, that in 
view of a rather small numerical value of F{(3) (equal to 
0.03 at the maximum), this enhancement factor may be 
not so big, despite a large value of the parameter rf/ 2 ql. 



III. SUMMARY 

We have presented an analytical study of the low- 
field magnetoresistance of a 2D electron gas subject to 
a weak one-dimensional modulation. We assumed that 
the disorder scattering is of small-angle nature due to the 
long-range character of the random potential. This cor- 
responds to the experimental situation in high-mobility 
semiconductor heterostructures, where the smoothness of 
the impurity potential is controlled by the large spacer 
separating the doping layer from the 2D electron gas. 

We have demonstrated that a strong magnetoresis- 
tance with a pronounced maximum at B ~ 0.4i? c (with 
B c defined by Eq. (||)) exists for sufficiently strong mod- 
ulation, rf^ 2 ql ^> 1. (In fact, due to smallness of a nu- 
merical coefficient, a rather large value of this parameter 
is required, rf/ 2 ql > 50.) The amplitude of Ap scales 
with the modulation strength r\ in this regime as 77 7 / 2 , 
see Eq. (|3^). At zero magnetic field the grating- induced 
correction to the resistivity is small and scales as 77 1 / 2 , see 
Eq. (p8|). At B » B c the magnetoresistivity is described 
by the earlier theory using the perturbative expansion in 
77 (Ref. U for white-noise random potential and Ref. ^ for 
smooth disorder). 

We have presented a detailed discussion of the limits 
of validity of the theory and, in particular, its matching 
with the earlier results of the 77-expansionQ Specifically, 
for a sufficiently weak modulation, rj 3 ' 2 ql <C 1, the re- 
sult of the 77-expansion (predicting Ap oc ij 2 ) is valid in 
the whole range of magnetic fields including the low-field 
region B < B c , implying disappearance of the strong low- 
field magnetoresistance. Therefore, the rf/ 2 behavior of 
Ap at B < B c , Eq. fl36|), does not imply a true non- 
analyticity of Ap(r?) at 77 — ► but rather restricts the 
applicability of the 77-expansion at weak B to the region 



G 



of sufficiently small modulation amplitudes, rf^ql <C i. 

Finally, to illustrate the magnitude of the effect, 
let us calculate A.p/p for typical experimental param- 
eters. Specifically, we use the parameters of a recent 
experimental electron density n — 2.84 x iO 11 cm~ 2 , 
modulation period a = 120nm, transport mean free path 
I = 19 fim. While the emphasis intl was put on a novel 
type of quantum magnetooscillations, a pronounced mag- 
netoresistance was observed in low magnetic fields, with 
a maximum at B m = 0.145 T. Using B m ~ 0AB c , 
we infer the modulation amplitude r\ = 0.157, so that 
rf' 2 ql ~ 62 is sufficiently large. With these values of 
parameters we find from (pq) the value of the magne- 
toresistivity at maximum, (Ap^/po)max — 6.8, in good 
agreement with the experimentally observed magnitude 
of the effect, p xx (B m ) / p xx (0) ~ 8. As concerns the resis- 
tivity correction at zero magnetic field, Eq. (|38l ) , we find 
Apxx(B — 0)/po — 0.49. Therefore, despite the para- 
metric smallness of this correction at 77 <C 1, its value 
in real experiments can be quite appreciable, due to its 
square-root dependence on the modulation strength rj. 
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APPENDIX A: RESISTIVITY CORRECTION AT 

B = 

In zero magnetic field and in the absence of disorder 
equations (g), (||) trivially decouple. We denote the ener- 
gies corresponding to the motion along x and y as E x and 
E y , respectively (E x + E y = Ep). The motion along x is 
unbound for E x > r]Ep: in the opposite case, E x < r]Ep, 
the particle is in a channeled orbit. The trajectories 
can be conveniently labeled by an angle ip such that 
E x — Ep sin 2 ip, E y — Ep cos 2 ip. Transforming the uni- 
form distribution on the Fermi surface (%,<fi) — 1 /2ira 
to the variables (x, ip) and integrating over x, we find the 
equilibrium distribution in the i/'-space, 



where 
r±(ip) = 

R{0,il>)-. 



R^^iP) , 



2vr 

dip 



1 



sin 2 ip 



cosz? 



1/2 



(Al) 

(A2) 
(A3) 



and ■& = qx. Eq. ( |A2] ) is valid for non-channeled orbits, 
sin 2 ip > rj (consideration of which will be sufficient for 
us to determine the conductivity, as explained below). 



The next step is the derivation of a kinetic equation 
characterizing the relaxation of n^(ip) toward equilib- 
rium due to disorder. The starting point is the Liouville- 
Boltzmann equation for n^,(x, (f>), 



1 

dtn<f, = -(pd^riA - idx-n^ + -5, 

T 



(A4) 



where the last term describes the diffusion in the velocity 
space induced by the small-angle scattering (in the ab- 
sence of this term ip would be an integral of m otion, im- 
plying dtriip — 0). Now we transform Eq. (A4) from the 
variables (x, <fi) to (x, ip) and use the fact that the oscilla- 
tions due to the motion in x-direction can be considered 
as a fast process as compared to the impurity s catte ring 
[the corresponding condition is specified in Sec. [II E| , see 
eq. (|^j)]. This allows us to average over x, which yields 
the following Fokker-Planck equation for n^{ip), 



(A5) 



In order to calculate the conductivity, we use the clas- 
sical Kubo formula, 



a xx = e v j dt 
'0 



dx 



2 7T 



a dx[ r 2n dcj/_ 
11 a J 2ti- 



x P(x, (j>; a/, </>'; t)v(x) sin <j>v(x') sin <p' , (A6) 

where v(x) — [2(Ep — V{x))/m] 1 / 2 is the particle velocity 
and P(x, (p; x' , cp'\ t) is the propagator in the phase space 
(x, 4>) (i.e. the probability density to move from a point 
{x 1 , (p 1 ) to a point (x, <p) in a time t). Transforming from 
the variables (x,<ft) to (r,ip), where r is the time along 
the trajectory, we get 



o , , .--.(- v j ill j dip 



dr- 



dip' 
27 



TWO 



dr' 



x v x (ip,T)v x (ip',T')vpsmip'P(T, ip;r', ip';t) , (A7) 

where T(ip) — — ° n ^ r_ (ip) is the period of oscillations 
in v x induced by the modulation. Averaging over the fast 



variable r, we reduce Eq. (A7) to the form 

dip' sin ip sin ip' 



e vvr 



dt I dip 



2ir r_ (ip) 



P(i/j,iP',t) 



(A8) 



The evolu tion kernel P(ijj,ip',t) satisfies the differential 
equation ( |A5[ ) with the initial condition P(ip,ip',0) = 
S(ip - ip'). Defining n(tp) = f™ dt J dtp'P(ip, ip',t) simp' 
and fi(ip) = n(ip)/r-(ip), we transform Eq. (A8) to the 
following form: 



■ > 1 dip , . 
(T. rj =e vv F j —smipn(ip) , 



-d^r + (ip)d^h(ip) = -sinip . 

T 



(A9) 
(A10) 
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The range of variation of the variable ip in (|A9|), (A1C) 
is restricted by the condition of non-channeled motion, 
| sin "0 1 > r] 1 / 2 . As soon as the particle comes into the 
region of channeled orbits, its velocity v x starts to oscil- 
late rapidly around zero, so that the contribution of such 
trajectories to the conductivity can be neglected (the cor- 
responding parameter is given by Eq. (]42|)). Therefore, 
Eq. ( A1C ) should be supplemented by the condition n — 
at the boundary of the non-chan neled region. In other 
words, the integral f(dip/2Tr) in (A9) is understood as 

In 1 



TT — 7] 



1/2 (dip/n) (we used rj -C 1), and the boundary con- 

3q. flAlcl ) reads n(ip = rj 1 / 2 , tt - rj 1 / 2 ) = 0. 
The function r + (ip) is equal to r + (ip) ~ 1 — jqV/ sm2 ^ 
for sin 2 ijj 77, reaching the value 2\/2/n at the bound- 
ary, sin 2 ijj = rj. In the limit 77 — ^ we have r + = 1 
and n(ip) = rsinip, yielding the Drude conductivity 
&xx = e 2 WpT/2 = (To. We now want to calculate the 
leading correction. It will be shown below to be of order 
rj 1 / 2 , so that we will neglect all contributions of higher 
orders. 

Let us denote by D the differential operator entering 
Eq. ( A10 ), D = -dipr + (->p)d^ on the interval [r] 1 / 2 ,^ — 
rj 1 / 2 ] with zero boundary conditions, so that Eq. (AG) 
can be rewritten as 



K--H 1 ' 2 fill; 

— sixiTpD- 1 shop . (All) 

1/2 7T 



Let us further consider a complete set of normalized func- 
tions on this interval (eigenfunctions of c* 2 ,) 



/»W0 



TT - 277 1 / 2 



1/2 . n(V> - rj 1 ' 2 ) 

1 - (2/71)^/2 ; 

n= 1,2,... . 



(A12) 



It is easy to see that to the first order in r] 1 / 2 we can 
insert the projector |/i)(/i| in Eq. (All), 



a xx = e 2 iyv 2 F r-(sm^\f 1 ) 2 (f 1 \D- 1 \f 1 ) 

7T 



(A13) 



where 



(/Is) 



di/>fW)g{il>) 



With the same accuracy we have 

(/ 1 |^- 1 |/i) = (M^I/i>- 1 



(A14) 



Expanding everything to the first order in 77 1 / 2 , we find 



(f 1 \D\f 1 ) = l + -(l-C 1 )r ! 1 / 2 + 0(r i ) 

7T 
1 /2 

(MsinV) = g) +0(rj), 



(A15) 
(A16) 



where 

./1 . n 2ir 



1-1- 



cosi? 



1/2' 



0.0241 



Substituting flAl4|) , ( |A15[ ), (|A16| ) in ( |A13[ ), we finally get 
a xx ^ao[l-C 2 V 1/2 ] , (A17) 

with a numerical coefficient C2 = -(1 — Ci) ~ 1.24. Since 
at zero magnetic field p^a; = cr";,. 1 , we arrive at the result 
( |38| ) for the grating-induced correction to resistivity. 
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